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Abstract 

An algorithm to compute the six distances between particles of 
a planar Four-Body central configuration is presented according to 
the following schema. An orthocentric tetrahedron is computed as a 
function of given masses. Each mass is placed at the corresponding 
vertex of the tetrahedron. The center of mass (and orthocenter) of the 
tetrahedron is at the origin of coordinates. The tetrahedron is orien- 
tated in a particular position function of the masses: with one of the 
particles placed on axis 3. The tetrahedron is rotated by two angles 
(to be tuned variables) around the center of mass until a direction 
orthogonal to the plane of configuration coincides with axis 3. The 
four coordinates of the vertices of the tetrahedron along this direction 
are identified with the weighted directed areas of the central config- 
uration. The central configuration corresponding to these weighted 
directed areas is computed giving rise to four masses and correspond- 
ing distances. The given masses are compared with the computed 
ones. The two angles of the rotation are tuned until the given masses 
coincide with the computed. The corresponding distances of this last 
computation determine the central configuration. The case with two 
equal masses also is considered. 

Keywords: Four-Body Problem. Planar central configuration. 

PACS 45.50.Pk Celestial mechanics 95.10.Ce Celestial mechan- 
ics(including n-body problems) 
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1 Introduction 



The coordinate system used in this paper has been presented in ref- 
erence pQ. The algorithm to compute planar central configurations 
starting from the weighted areas was published with several applica- 
tions in [2]. 

To begin I make a review of the new four body coordinates that 
mimics an important portion of [I], expanded at important points and 
reduced in several comments. 

The masses of the four bodies mi, 7712, 1TI3 and 7714 are positive, 
generally different, but the values could be repeated. 

We transform from the inertial referential, to the frame of principal 
axes of inertia by means of a three dimensional rotation G parame- 
terized by three independent coordinates. 

In addition to this rotation three more coordinates are introduced, 
as scale factors R\, R2, -R3, which are three directed distances closely 
related to the three principal inertia moments through 

h = li{R 2 2 +Rl), I 2 = »{Rl + R 2 1 ), and h = l i{R\ + R 2 2 ), (1) 

where /j is the mass 



_ 3/ mi m 2 m 3 m 4 
V mi + ni2 + 77J3 + 7744 

With the first rotation and the change of scale the resulting four- 
body configuration has a moment of inertia tensor with the three 
principal moments of inertia equal. The second rotation G' does not 
change this property. 

The cartesian coordinates of the four particles, with the center of 
gravity at the origin, written in terms of the new coordinates are 

/ a\ a>2 03 04 
h b 2 63 h 
\ ci c 2 c 3 c 4 
(3) 

where G and G' are two rotation matrices, each one a function of 
three independent coordinates such as the Euler angles, and where 
the column elements of the constant matrix 

ai 02 03 04 
E = I 61 b 2 h 64 I ( I ) 

Cl c 2 c 3 c 4 
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are the coordinates of the four vertices of a rigid orthocentric tetrahe- 
dron, with the center of mass at the origin of coordinates, namely: 



dimi + a2?7i2 + 037713 + (247714 = . 
bimi + b 2 m 2 + 637773 + 647714 = , 
c\m\ + c 2 m 2 + c 3 m 3 + C47714 = . 



(5) 



We introduce the following notation for the matrix 



M = 



/ mi 





777 2 









7773 













7774 / 



\ 



(6) 



V 



An equivalent condition in order to have three equal inertia mo- 
ments for the rigid tetrahedron is expressed as 



An orthocentric tetrahedron has the property that the perpendic- 
ular lines to the faces trough the four vertices intersect at the same 
point. Orthocentric tetrahedra were considered by Lagrange in 1773 
[3] . Other old references on orthocentric tetrahedra are found in a pa- 
per by Court [1], he calls them orthocentric and orthogonal to these 
tetrahedra. Placing the four masses at the corresponding vertices, the 
intersection point is actually the center of mass of the four masses, 
and the moment of inertia tensor of the four particles has the same 
principal value in any direction. Equation ([7]) imply that the inertia 
tensor of the rigid tetrahedron is proportional by a factor 2/_t to the 
unit matrix. 

To show these properties we consider for a while a fourth coordi- 
nate for the four vertex of the tetrahedron 



where we use the notation 777 = rrt\ + 7712 + 7773 + 7774 for the total mass 
of the system. Then using properties ([5]) , (JT]) , and © we write them 



/ 1 \ 
EME T = /j, 10 

V° 1/ 



(7) 
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in terms of r 



— in the form 



1 



( 0,1 a<2 a 3 04 \ 

61 b 2 63 b A 

C\ c 2 c 3 c 4 

\ r r r r ) 



M 



( ai 61 ci r ^ 

a 2 62 c 2 r 

a 3 63 c 3 r 

y a 4 64 C4 r y 



/ 1 \ 

"0100 

10 

V 1 j 



(9) 

Since the inverse matrix from the left is equal to the inverse from the 
right, this equation transforms into 



/ dl 


h 


Cl 


r 




/ 


ai 


02 


03 


04 








( 


mi 











\ 




b 2 


C2 


r 






h 


b 2 


63 


64 













JL. 

m 2 










«3 




C3 


r 






Cl 


c 2 


C3 


c 4 
















JL. 

m 3 







\ 04 


64 


C4 


r 


/ 


V 


r 


r 


r 


r 


J 






\ 











JL. 

m,4 


J 
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This matrix equation is equal to its transposed; therefore it just has 
ten independent equations. Four of them are 



+ <? = M [J__1] (j = 1,2,3,4), (11) 



mo- 



rn 



for the square of the distance from the orthocenter to the vertex j. 
The other six are 



ciiCij + bibj + CiCj = -— (i^j), 



m 



(12) 



the internal product of two directed vertices is the same for the six 
different possible choices of the pair From these basic equations 
it is easy to show that the position vector of one vertex is orthogonal 
to the three vectors between two vertices of the oposite face (to the 
position vector of the vertex.) 

ai(a,j - a k ) + bi(bj - b k ) + Ci(cj - c k ) = (i,j,k different.) (13) 

In addition, the square of the distance between two vertices is given 
by 

- aj) 2 + (bi - bj) 2 + - cj) 2 



TO; 



TOo' 



This last is the condition to have a moment of inertia tensor with 
the same three principal moments of inertia. The six edges of the 
tetrahedron should be equal (for an arbitrary (i) to the square root of 
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the right-hand side of this equation. The volume of this tetrahedron 
is equal to 1/6 if fx is selected as above. 

There are other remarkable geometrical properties of an orthocen- 
tric tetrahedron. The center of mass of each face is at the orthocenter 
where the three altitudes of the face intersect. This point is on the 
same straight line between the opposite vertex and the center of mass. 
In addition to the orthogonality of the three sets of two opposite edges 
of the tetrahedron, the two orthogonal edges are also orthogonal to 
the line joining the center of mass of the two edges. 

There are important coordinate systems to fix the origin for mea- 
suring the G' rotation; from these I prefer to choose in this paper 
particle with mass mi along coordinate axis 3, the other three in a 
parallel plane to the coordinate plane of 1 and 2, which does not in- 
clude the first particle; the particle with mass m,2 on an orthogonal 
coordinate plane that includes the first particle and the center of mass, 
and the other two particles on a line that is parallel to the coordinate 
axis 1 and perpendicular to the coordinate plane of the first two par- 
ticles. Particle 1 has the coordinates 



(«L.h.n)= \0,0,\ Km — I ( ir >) 



Particle 2 has the coordinates 



(«.*,■<»)-( 0, J " ( 7 + mi \ , -Jr^V ■ (16) 
V y m2{m — mi) y (m — m\)m J 

Particle 3 has the coordinates 

(a 3 ,& 3 ,c 3 ) = 



/iTIl4 / /i77l2 / fimi 



\ y m 3( m 3 + m ±) ' y ( m 3 + mi)(m — mi) ' y (m — mi)m J 

(17) 

Particle 4 has the coordinates 

(04, &4, C4) = 



/J7713 / fim2 / [imi 



777,4(777.3 + 777,4)' y (7773 + 7774) (m — mi ) ' y (777 — 7771)777/ 

(18) 
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This rigid tetrahedron is the generalization of the rigid triangle of 
the Three-Body problem with the center of mass at the orthocenter 
discussed previously in [5]. The same triangle was used with different 
purposes by C. Simo [6]. 

I assume for simplicity that the potential energy is given by the 
Newton potential (the gravitational constant is equal to 1) 



V 



i<j 



(19) 



although our results may be generalized for any potential with a given 
power law of the relative distances between particles ry . It follows the 
relation between the inter-particle distance and the new coordinates. 
The relative position between particles i and j is 



/ 










Vj 






V 


Z 3 


-Zi ) 






( 



\ C 3 



hi 

Ci 



(20) 



The square of this vector is not a function of the first rotation G, but 
just of the scale matrix and the second rotation matrix 

/ a* — a 



(21) 



where A is the symmetric matrix 

A = A 12 
V A 13 




G' 












G' J 



(22) 



The six distances are thus functions of six components of matrix A or 
equivalently, are functions of the six independent coordinates in the 
scales Ri, and the rotation G'. 

We also compute the kinetic energy as a function of the new coor- 
dinates, which is given by 



2 



E^ 2 



4(R 2 R 3 u} 1 n 1 + R 3 Rmj 2 ^2 + RiR 2 ^3) + 



=i 



00 



( R 2 + i?§ 





Rl + R\ 




0J + 



6 



( R\ 



V 



o 
o 







R 2 ' » 2 







R\ + R 2 2 



(23) 



where w = (wi,W2i w 3) is the angular velocity vector of the first ro- 
tation G, and £1 = (Q,i,£l 2 ,tl 3 ) is the corresponding angular velocity 
vector of the second rotation G'. 



2 Equations of Motion 

The equations of motion follow from the Lagrange equations derived 
from the Lagrangian T — V as presented in any standard text on 
Mechanics [TJ, 0. 

However, the three coordinates related to the first rotation pro- 
duces Lagrange equations that imply, when the potential energy is a 
function only of the distances, conservation of the angular momentum 
vector in the inertial system 

GL, (24) 

where L is the angular momentum in the principal moments of inertia 
frame 

L - — 

duo 

( {R 2 2 + Rl)u x \ ( R 2 R 3 n 1 

(Rl + R\)u2 - 2/J R 3 Ri$l 2 I . (25) 



^ (Rf + Rl)u> 3 J \ R ± R 2 n 
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This conservation leads to three first order equations forming, for 
this four-body problem a generalization of the Euler equations (valid 
for the rotation of a rigid body), namely 



d 
~dl 



I n(R% + RI)uji - 2fiR 2 R 3 n 1 
fj,(Rl + R\)lo 2 - 2fiR 3 R 1 Q 2 
\ niRl + R%)u) 3 - 2/i J R 1 ^ 2 Q 3 



/ n(R\ - Ri)uj 2 U3 + 2fiRi(R 2 uj 2 n 3 - R 3 oj 3 ^2) 

- B%)u 3 ui + 2fiR 2 (R3U} 3 ni - Riuii^s) | . (26) 
\ n(Rl - Ri)u>iu) 2 + 2fj,R 3 (R 1 u} 1 n 2 - i? 2 a>2^i) 

The so called elimination of the nodes in the Three-Body problem 
, has a similar representation in this coordinates for the Four-Body 
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problem by means of the equation that equals the angular momentum 
vector in the principal moments of inertia frame to the rotation of a 
constant vector, which may be written in terms of two Euler angles 

/ (Rl + B$) Ul \ ( R 2 R 3 n 1 \ ( \ 

A* {Rl + Ri)^ 2 - 2/x R 3 RiSl 2 = £G T , (27) 
\ (Rl + R 2 2 )co 3 J \ R!R 2 n 3 J \ 1 J 

where I is the magnitude of the conserved angular momentum. 

The Lagrangian equations of motion for the three scale coordinates 

are 

d 2 „ _ r „ „ „ „ , „ . , , „o „o, dV 



Vj^Ri + 2fi[R 2 oj 3 n 3 + R 3 u 2 n 2 ] + + w| + ^ + ng) - . 

(28) 

Ai^2 J R2 + 2 / ti[-R 3 Wl^l + J Rl^3^3]+^2(wi+W 2 + ^i + ^?) = —Qjj-: 

(29) 

and 

d 2 <9F 
Ai— ^i?3 + 2/x[i?io;2^2 + i?2Wifii] + V-R^l + ul + ti{ + = . 

(30) 

The three equations of motion for the three coordinates associated 
with the second rotation G' are written as an Euler equation similar 
to the one found for the first rotation, although the internal angular 
momentum is not conserved because of the presence of an internal 
torque 



d_ 
~dt 



( n{Rl + Rj)^ - 2fiR 2 R 3 uj 1 
fi(Rl + Rj)n 2 - 2 t iR 3 R 1 co 2 
\ n{R\ + R\)VL 3 - 2nRiR 2 u 3 




( n(R% - Rj)n 2 n 3 - 2nR 1 (R 2 u 2 n 3 - R 3 u 3 n 2 ) \ 

fi(R 2 - J?|)0 3 ni - 2nR 2 (R 3 u> 3 £li - R^n 



3) 



(31) 

y n(Rl - B%)£l£h - 2(iR 3 (R 1 u 1 tt 2 - ifcc^fii) j 

where K±, K 2 , K 3 are the components of the internal torque K which 
is expressed in terms of the derivatives of the potential energy with 
respect to the three independent coordinates qj in the rotation G' 
and the three vectors Cj that appear in the expression of the angular 
velocity in terms of the same coordinates 

3 

n = (32) 
3=1 
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where the vectors Cj are generally functions of the coordinates qj . 
The internal torque is determined by the equations 



K c, 



dV 



(33) 



There is one more constant of motion, namely the total energy 

3 



E = T+K = V+- 
2 



~ 4(^2-^3^1^! + R 3 Rl^2 + RlR 2 0J3^3) + 



U) 



T 



r 2 

i=l 

(Rl + Rl 

Rl + Ri 

Rl + Rl 

Rl + Rl 

Rl + Ri 

Rl + Rl 



UJ + 



n 



(34) 



3 The plane problem 

The case with the four particles in a constant plane is an important 
and old subject [10J. Our coordinates are now adapted to that case. 
The third components of the cartesian coordinates of the four particles 
are zero. The modification of our coordinates (4) for this case is given 
by two changes: the first one is a rotation by just one angle in the 
plane of motion; and the second, the scale associated with the third 
coordinate is set equal to zero, namely 



/ X\ X2 X3 X4 
Vi V2 V3 Vi 







^ cos tp — sin ip 
sin tp cos ip 
v 1 



This equation simplifies to 




/ ai d2 03 04 
h b 2 b 3 6 4 
\ ci c 2 c 3 c 4 



(35) 



Xl X2 X3 X4 

yi V2 V3 Vi 
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/ cosV> -simp \ ( Ri \ ,t I ? X ? 2 ? 3 ? 4 1 

V 7 V 7 \ Cl C 2 C 3 C 4 / 

(36) 

in terms of six degrees of freedom. 

We need three independent coordinates (for example three Eu- 
ler angles) in G' for the two independent vectors in four dimensions 
expressed in the base of the three constant vectors a, b, and c, or- 
thogonal to the mass vector. 

In order to formulate in a mathematical language the conditions for 
a plane solution we use the equations introduced by Dziobek [10] , that 
have been promoted by many years by A. Albouy and coworkers (see 
|llj and references therein,) which consists in using the four directed 
areas of the triangles formed by the particles. 

The four (double) directed areas are written in terms of the carte- 
sian coordinates as 





1 


1 


1 




1 


1 


1 


Si = 








, s 2 = 


Xx 


X i 


£3 




2/2 


2/3 


2/4 




2/1 


2/4 


2/3 




1 


1 


1 




1 


1 


1 


S 3 = 




X2 




, S4 = 


Xx 




X2 




yi 


2/2 


2/4 




2/1 


2/3 


2/2 



which are the four signed 3x3 minors formed from the matrix 

/ 1 1 1 1 \ 

xx X2 x 3 x 4 . (38) 
\ 2/1 2/2 2/3 2/4 / 

From now on I will refer to these areas simply as the directed areas. 

The addition to the previous matrix of a row equal to any of its 
three rows produces a square matrix with determinant zero, that im- 
plies that the necessary and sufficient conditions to have a constant 
plane tetrahedron are 

4 

E^ = ' ( 39 ) 
i=l 

and 

4 4 

E^ = o> E^ = ' ( 4 °) 
i=i i=i 
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The two last equations are summarized by the zero vector condition 



i=i 



(41) 



An expression for the three directed areas in terms of the previous 
coordinates follows 



s 3 



CME T G 




(42) 



where C is a constant with units of area over mass. Equation f|39|) is 
satisfied from this expression of the directed areas, upon using that the 
four-vectors in matrix E are orthogonal to the mass vector. With the 
substitution of equations (f36|) and (|4"2|) in equations ([4*0]) one obtains 
an identity, independent of coordinates R\, R 2 , tp, and of the rotation 
angles of G'. 

I denote here 9 and cj> the spherical coordinates determining this 
vector. 



(43) 



Given the four masses, the four directed areas of the four particles 
are functions of this unit vector direction only, up to a multiplicative 
constant C depending on the choice of physical units. Other form 
equivalent to the plane condition (j41[) is also published in reference 









' sin 9 cos <p \ 






- 


sin 9 sin (j) 








1 cos 9 J 



0. 

In the plane case the angular momentum has a constant direction 
orthogonal to the plane and of magnitude 

p, = ^ 

dip 

The kinetic energy becomes 
" 2 

£ Ri - 4(R 1 R 2 ipn 3 ) + i> 2 (Rl + Rl) + 



n[ip(Rj + Rl)-2R 1 R 2 n 3 



(44) 



T 



2 



/ Rl 






(45) 
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Now I substitute polar coordinates for the R\ and R2 coordinates 

R\ = R cos a , i?2 = R sin a . (46) 

Writing the kinetic energy in terms of the angular momentum constant 
of motion instead of the ip velocity lead us to 



R 2 + R 2 (a 2 + J2| cos 2 (2a) + ft 2 sin 2 a + n| cos 2 a 



p2 
1 



The energy conservation is thus expressed as 



E 



i? 2 + i? 2 d 2 + ft 2 . cos 2 (2a) + n; sin 2 a + fi| cos 2 a 



2^i? 2 ' 
(47) 



P 2 



where V represents the potential energy. 



2fiR 2 
(48) 



4 Central configurations 

The non collinear planar central configurations are characterized in our 
coordinates by constant values of the G' matrix and of the coordinate 
a associated to the constant value of the ratio R\/R2- For these cases 
the angular velocity vector Q, is the null vector, the angular velocity 
d is also zero and the conservation equations of moment and energy, 
([31]) and ([15]) respectively, become 

= ^R 2 . (49) 

and 

E =i k2+ S? +v - <5o) 

Note that the potential function V becomes equal to a constant di- 
vided by coordinate R. 

These equations are identical to similar equations obtained for the 
Euler and Lagrange central configurations of the Three-Body problem 
|12j . They are formally the same as the equations for the conies in the 
Two-Body problem in terms of the radius R and the true anomaly ij). 

The constant values of the G' matrix and angle a referred to above 
are not arbitrary but they are determined by three independent quan- 
tities as discussed in what follows. 
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I begin with the approach by Dziobek [10] . The Four-Body cen- 
tral configurations are determined as critical points of the potential 
energy with a fixed total inertia moment. The planar solutions with 
zero volume but finite area are obtained taking in account that the 
variational equation is modified adding the restriction of planar mo- 
tion. This condition is obtained by Dziobek [10J from the derivative 
of the Cayley-Menger determinant with respect to rfj that he found 
to be proportional to the product of the directed areas SiSj. 

It follows that the solution is given in terms of Lagrange's multi- 
pliers A and a 



_-3 
r jk 



a + XAjAf. . 



(51) 



where Aj = Sj/rrij are weighted areas, quotient of the directed area 
divided by the corresponding mass. This equation was presented by 
Dziobek [10j . A proof was published by Moeckel [15] : and using a 
different approach to the same problem, deduced by Albouy [IT] . A 
new proof of the equation was obtained in a different approach by 
Piha and Lonngi [2J. 

It follows from (|42p that in a planar solution the weighted directed 
areas are expressed as 



(Ax \ 

A 2 

^3 

V A, J 



CE T G' 



/ 



sin a cos ( 
sin 9 sin ( 
cos 9 



(52) 



The weighted directed areas are up to a normalization factor equal 
to the third rotated coordinate of the rigid tetrahedra. In terms of 
vectors (15-18) and angles 9 and <f> this equation is expressed as 



Aj = C {aj sin 9 cos <j) + bj sin 9 sin (j) + Cj cos 9) . 



(53) 



The explicit expressions choosing C = \J {m — mi)/ n are 

m — ml 
Ai = — -—=■ cos 9 , 

Jmim 



A 1711 a , m 3 + m 4, ■ a ■ 
A2 = — \ cos 9 + \ sin 9 sm< 

V m y 7712 



a a I ^ ■ a ■ a, / m ±( m -mi) 

A3 = —\ cos 9— , / sin sin <p+\ z r sint^ cos 1 

V 777 V m 3 + m 4 V m 3\ m 3 + m <±) 
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/mf „ / m2 •/,„•- / m 3 (m-mi) . 

A4 = —\ — cos u—\ sin t) smd)—\ r- sin cos <p . 

V m V m 3 + m 4 V m 4\ m 3 + 777 -4j 

(54) 

Since the lengths and masses are defined up to arbitrary units, we 
assume [2], with no loss of generality, that the parameter a equals 
unity 

r^ = l + XAjA k (j^k). (55) 

This implies, according to [2J that the unit of distance always obeys 
the restriction 

E rrijm k /r jk 
a = ^ - = 1. (56) 

j>k 

In the paper by Piha and Lonngi [2 J the directed weighted areas are 
known as four given constants. The previous equation then gives the 
distances as functions of the unknown parameter A. Through them, 
the areas of the four triangles become functions of A, that should obey 
the necessary restrictions (139|) . (|40j) . to verify that one has a planar 
solution. This restriction allows in many cases to determine the value 
of A and hence the values of the six distances and the four masses. 
This is an implicit way to deduce planar central configurations with 
four masses. 

In this paper we assume the four masses are known from the begin- 
ning, the four weighted areas are then determined by expressions (|54p 
in terms of the two tuning variables 9 and <j). Given particular values 
of these two angles, become determined the four constants Aj, until a 
multiplicative factor, which in turn produces a computed central con- 
figuration with given distances and masses. These computed masses 
are not in general equal (or proportional) to the starting values used to 
compute the orthocentric tetrahedron. The two angles are then tuned 
until a numerical coincidence is produced between the given and the 
computed masses. The distances between particles, computed for this 
last central configuration, correspond to the given masses. 

Choosing the values of the masses as mi = 10, 771,2 = 13, 7773 = 
15, 7714 = 17, I have been able to compute seven central configurations. 
A concave with mi at the interior for A = —2.32656490060845 

6 = 0.139240050165164 <f> = 4.8453912490189 

r 3 i = 0.749352173668766 r 42 = 1.28360719430431 

r 4 i = 0.730065912777101 r 23 = 1.22971324110202 

= 0.639643905964532 r 43 = 1.10420199339901 
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Figure 1: Spherographic projection of the hemisphere of the two angles mo- 
tion of the orthocentric tetrahedron. The great circles represent the positions 
where three particles are collinear. The four spherical triangles are concave 
open sets labeled by the particle at the interior of the triangle. The spherical 
rectangular open sets correspond to convex configurations with the same or- 
der that the neighbor triangles. The isolated points are at the angles where 
a central configuration has been computed. The values of the masses are 
mi = 10,m2 = 13, = 15, = 17. 
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A concave with m 2 at the interior for A = -1.82536262997443 



6 = 1.03401714345126 <f> = 

r 23 = 0.804939324544205 r 4 i 

r 42 = 0.798422993763499 r 31 

ris = 0.640524087661698 r 43 



4.70428233126543 
= 1.36218205283364 
= 1.32975926553097 
= 1.07615875773082 



A concave with m 3 at the interior for A = —1.6355867896085 



6 = 1.00916705438091 <f> = 

r 23 = 0.819159125672848 r 4 i 

r 43 = 0.807449041118773 r 12 

r 31 = 0.666344917442357 r 42 



2.69053445366272 
= 1.41307775730063 
= 1.34385908526924 
= 1.08738922587974 



A concave with m 4 at the interior for A = -1.55004225977105 



6 = 1.33726059387504 = 

r 4i = 0.80384936609868 r 23 

r 43 = 0.768450306675012 r i2 

r 42 = 0.730824120911527 r 3 i 



= 0.589297944199426 
= 1.44896029979243 
= 1.27362211521548 
= 1.18407793737811 



A convex with r 4i at the diagonal for A = -0.59067068058041 



6 = 0.833903746254753 
r 4i = 1.26246852646001 
ri2 = 0.834421340108319 
r 31 = 0.855040198605932 



= 3.74619198490858 
r 23 = 1.21353489669217 
r 43 = 0.914579880279118 
r 42 = 0.900505139131835 



A convex with r 42 at the diagonal for A = —0.571602583650311 



6 = 0.849784880696665 
r 42 = 1.1712284103105 
ri2 = 0.837603049326044 
r 4i = 0.874842967906726 



= 5.66013081680728 
r 3 i = 1.31084659340392 
r 43 = 0.916446768964837 
r 23 = 0.888555880648135 



A convex with r 43 at the diagonal for A = -0.532931485997706 



6 = 0.861931053448714 = 

r 43 = 1.13115400784062 r 12 

r 3 i = 0.863912312413253 r 42 

r 4i = 0.880188793100079 r 23 



1.64118840218233 
= 1.36414262641605 
= 0.907050181187521 
= 0.893570787310162 
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5 The case with two equal masses 



When two masses are equal, the other two generally different, the 
orthocentric tetrahedron become symmetric with respect to a sym- 
metry plane containing the two (possible) different masses and the 
center of mass. When the tetrahedron is projected orthogonal to a 
direction contained in the symmetry plane the configuration becomes 
with the kite symmetry. Central configurations with two equal masses 
and kite symmetry are allowed if and only if the plane of symmetry 
is projected in a line. Such particular case is considered rotating the 
symmetry plane by a single tuning angle until the given masses coin- 
cide with the computed ones. It is only relevant half a circle of this 
angle. Three cases are possible separated in the semicircle by three 
contiguous different open sectors. Two concave sectors and one convex 
sector. 

One illustrate this case computing central configurations for the 
masses mi = 8, m 2 = 10, m 3 = m 4 = 9. Two different symmetric 
concave central configurations were found with the largest mass in 
the interior for A = -1.938380387574535 

= 0.053846821928071 </> = vr/2 

r 12 = 0.771017850970647 r 43 = 1.37174365887618 

r 31 = 0.737687261978754 = r 41 

r 23 = 1.24798073364597 = r 42 

and for A = -1.9051833373390781 

6 = 0.131677055692469 <f> = 3^/2 

ri2 = 0.704334104242389 r 43 = 1.16065690193433 

r 31 = 0.777736225967895 = r 4 i 

r 23 = 1.35289892174911 = r 42 

and one convex for A = -0.6629645165685381 

6 = 1.00225261834462 <f> = tt/2 

r 43 = 1.2388066072406 r 12 = 1.24313752444539 

r 3i = 0.891030948126205 = r 4 i 

r 23 = 0.864175379988576 = r 42 . 

No concave solution was found for the lowest mass mi = 8 in the 
interior. 

Non symmetric central configurations do exist for the case of two 
equal masses. For the given masses mi = 8,m 2 = 10, m 3 = m 4 = 9 
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one has the non symmetric convex central configuration obtained 
by searching with two angles as tuning variables with r 23 at the 
diagonal, and the two equal masses not on a single diagonal, for 
A = -0.65459120484497 

9 = 0.991809150592805 <p = 3.71189655731751 

r 23 = 1.28383284893302 r 4i = 1.20043491498806 

r 42 = 0.865256456835286 r 31 = 0.891936018662304 

ri2 = 0.878638767853787 r 43 = 0.879635445676687 

and the the non symmetric concave central configuration with one of 
the equal masses at the interior for A = —2.0782713344769 

6 = 1.2731702593358 <f> = 2.4920236726443 

r 23 = 0.681773343374296 r A1 = 1.13338465823949 

r 31 = 0.767118374257656 r 42 = 1.31028169747428 

r 43 = 0.776744110452106 r l2 = 1.35194750927076 
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